Damage spreading in the mode-coupling equations for glasses 
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We examine the problem of damage spreading in the off-equilibrium mode coupling equations. The 
study is done for the spherical p-spin model introduced by Crisanti, Horner and Sommers. For 
p > 2 we show the existence of a temperature transition To well above any relevant thermodynamic 
transition temperature. Above To the asymptotic damage decays to zero while below To it decays 
to a finite value independent of the initial damage. This transition is stable in the presence of 
asymmetry in the interactions. We discuss the physical origin of this peculiar phase transition 
which occurs as a consequence of the non-linear coupling between the damage and the two-time 
correlation functions. 



The theoretical understanding of the dynamical be- 
havior of glasses is a long outstanding problem in sta- 
tistical physics which has recently revealed new aspects 
related to the underlying mechanism responsible of the 
glass transition [ljj^]. While there are still some obscure 
points in the theory (i.e. the inclusion of finite time acti- 
vated process beyond the mean-field limit) a scenario has 
emerged which unifies the dynamical approach (mode- 
coupling theory) with the thermodynamic Adam-Gibbs- 
Di Marzio approach. The scenario for the dynamical 
behavior of glasses can be summarized in three differ- 
ent temperatures which separate three different regimes. 
In the high-temperature regime T > Td the system be- 
haves as a liquid and is very well described by the mode- 
coupling equations of Gotze in the equilibrium regime || . 
A crossover takes place at Td where there is a dynamical 
singularity and the correlation functions do not decay 
to zero in the infinite time limit (ergodicity breaking). 
This dynamical singularity is a genuine mean-field effect 
which turns out to be a crossover temperature when ac- 
tivated processes are taken into account. Below T d the 
relaxation time (or viscosity) starts to grow dramatically 
fast and seems to diverge at T s where the configurational 
entropy vanishes and a thermodynamic phase transition 
takes place. The glass transition T g (as denned where the 
viscosity is 10 13 Poise) lies between T s and Td and de- 
pends on the cooling rate. Hence T g does not correspond 
to a true dynamical singularity. Furthermore, there is a 
first order phase transition Tm where the liquid (if cooled 
sufficiently slow) crystallizes. 

The essentials of this scenario have been corroborated 
in the context of mean- field spin glasses, and in particular 
in those models with a one step replica symmetry break- 
ing transition While the first-order transition tem- 
perature Tm is absent in spin glasses (disorder prevents 
the existence of a crystal state) the other two transitions 
(T s ,Td) have been clearly identified. 

The purpose of this paper is the study of the dam- 



age spreading in models for structural glasses. Dam- 
age spreading is the study of the time propagation of 
a perturbation or damage in the initial condition of a 
system. The propagation of the initial damage is a dy- 
namical effect which has deserved considerably attention 
in the past (specially in the context of dynamical sys- 
tems, for instance, networks of boolean automata |J) 
because it allows to explore the structure of the phase 
space of the system. To investigate damage spreading 
a suitable Hamming distance in the space of configu- 
rations is defined. Then we consider two random ini- 
tial configurations {(Ji,Ti} with a given initial distance 
Do for two identical systems which evolve under iden- 
tical noise realizations and compute the distance D(t) 
as a function of time. Eventually one is interested in 
the asymptotic long-time behavior of the distance D(t), 
i.e -Doo = limt—,00 D(t). In general, three different 
regimes can be distinguished. A high-temperature regime 
T > Tq where D^ = independently of the initial dis- 
tance Dq. A intermediate regime T\ < T < To where 
Doo = (T) is not zero but independent of the initial 
distance. And finally a low-temperature regime T < T\ 
where D^ = D^iT, D ) depends on both temperature 
and initial distance. While it is widely believed that T\ 
corresponds to a thermodynamic phase transition it is 
not clear what the physical meaning of To is. 

In this work we address the problem of damage spread- 
ing in the off-equilibrium mode-coupling equations. We 
show the existence of the temperature To in glasses well 
above T d and T s in the high temperature phase. We 
show that this new transition is a consequence of the 
non-linear coupling between the damage and the corre- 
sponding two-times correlation function. This effect is 
an essential ingredient of the mode-coupling equations 
and should be generally valid even beyond the mean- 
field limit. We believe the appearance of this damage 
transition is a quite general result in glassy models (with 
and without disorder) where the scenario of Gotze for 
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mode-coupling transitions is valid. 

The simplest solvable model which yields the off- 
equilibrium mode coupling equations is the spherical p- 
spin glass model In this case, the configurations are 
described by N continuous spin variables {af, 1 < i < N} 
which satisfy the spherical global constraint Yli=i °f = 
N. The Langevin dynamics of the model is given by, 



3(7-1 

— = Fi{{a}) - (J,(Ti + r}i 
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where Fi is the force acting on the spin cr, due to the 
interaction with the rest of the spins, 
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and TL is a Hamiltonian. The Jf ' "' ' " are 
quenched random variables with zero mean and variance 
p\/(2N p ~ 1 ) which we take to be symmetric under per- 
mutation of the different super indices. The calculations 
presented here can be easily generalized to asymmetric 
couplings 0. Obviously, in this last case there is no en- 
ergy TL which drives the system to thermal equilibrium. 
The term /i in eq.(|^) is a Lagrange multiplier which en- 
sures that the spherical constraint is satisfied at all times 
and the noise rj satisfies the fluctuation-dissipation rela- 
tion (i]i(t)rij(s)) = 2T5(t — s)Sij where (...) denotes the 
noise average. 

We define the overlap between two configurations of 
the spins a, r by the relation Q = YliLi a i T i an d the 
Hamming distance between these two configurations as, 
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in such a way that identical configurations have zero 
distance and opposite configurations have maximal dis- 
tance. Now we consider two copies of the system {er^, r^} 
which evolve under the same noise (0) but with differ- 
ent initial conditions. Here we restrict to random initial 
configurations (i.e equilibrium configurations at infinite 
temperature) with initial overlap Q(0). Nevertheless, dif- 
ferent type of initial conditions can be considered. The 
different set of correlation functions which describe the 
dynamics of the system are given by, 
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where hf,hj are fields coupled to the spins <Ji,T; L re- 
spectively. In what follows we take the convention t > s. 
The previous correlation functions satisfy the boundary 
conditions, C(t,t) = l,R(s,t) = 0, lim t _,.( s )+ R(t, s) = 1 
while the two replica overlap Q(t, s) defines the equal 
time overlap Qd{t) = Q(t, t) which yields the Hamming 
distance at equal times or damage D(t) through the re- 
lation (||). Following standard functional methods |||]] 
it is possible to write a closed set of equations for the 
previous correlation functions, 
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duR(t,u)C(s,u)C p - 2 (t,u) (7) 
+ n(t)R{t,s) = S(t -s) + 
duR(t,u)R(u,s)C p - 2 (t,u) (8) 



■ n(t)Q{t,s) = | / duR{s,u)Q p -\t,u) + 
A Jo 
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duR(t,u)Q(u,s)C p ~ 2 (t,u) (9) 



while the Lagrange multiplier n(t) and the diagonal 
correlation function Qd(t) obey the equations, 



l*(t)=T + ^ f duR{t,u)C p - 1 {t,u) (10) 
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duR(t,u)Q(t,u)C p ~ 2 (t,u) (11) 



This set of equations is quite involved. For the correla- 
tion C and response functions R eqs. (|t|J^,|io1) several re- 
sults are known, in particular their behavior in the FDT 
regime (where time translational invariance is satisfied 
and the fluctuation-dissipation theorem is obeyed) as well 
as in the aging regime || . 

A first glance to equations (||),([ll]) reveals that the 
overlap Q(t,s) and its diagonal part Qdit) are coupled 
each other through the correlation C(t,s) and response 
function R(t,s). The trivial solution Q(t,s) = C(t,s) 
and Qd{t) = 1 corresponds to the case where the ini- 
tial conditions are the same, Qd(0) = 1 and the distance 
D(t) = for all times. This solution (hereafter we will 
denote it by HT) is asymptotically reached by the dy- 
namics for high enough temperatures. The typical time 
needed to reach that solution grows if temperature de- 
creases. At a given temperature (which we identify with 
To) there is an instability in the dynamical equations 
© i (0) an d the asymptotic solution differs from the HT 
one. We did not succeed in finding an explicit expression 
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for To but we have been able to show its existence and 
find lower and upper bounds for its value. 

To show the existence of To we focus in the high tem- 
perature FDT regime (t — s)/t << 1 with t, s both large 



and TR(t, s) = TR(t - s) 



22j*r«i. Writing Q(t,a) = 
Qd(s)Q(t,s), using the inequalities Q(t,s) < C(t — s), 

9Q d (t) 



— Q- t — > and inserting these results in 
to get the following inequality, 
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A trivial solution which always satisfies that inequality 
is Q d (t) = 1. It is not difficult to check that the previous 
inequality yields a lower bound for the temperature at 
which there is an instability in the condition ( Il2| ) , leading 
to 



To > 
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On the other hand, a linear stability analysis of the equa- 
tions (|),(|l|) around the HT solution, Q d (t) = 1- ef(t), 
Q{t, s) = C(t -s)- eg(t, s) where /(0) = 1, g(t, t) = f(t) 
yields for equation (0) in the large t limit, 



dg(t, u) 
du 



(14) 



Finally, the inequality d9 Q^ > together with ( |l3| 
yields, 



- 1 < Tn < 



(15) 



As said before, it is very difficult to find an explicit 
expression for To. The reason is that both Qd(t) and 
Q(t, s) (or equivalently, f(t) and g(t, s) in the linear sta- 
bility analysis) are not related by any TZ?T-like relation 
in the long time limit. Consequently, the analysis of the 
dynamic instability turns out to be more difficult. 

Note that for the particular case p = 2 the inequality 
( |l5| ) yields To < 1. Taking into account that (]l5|) was 
derived under the assumption To > T s = 1 that yields 
To = 1 . The simpler case p = 2 has been already consid- 
ered by Stariolo |l(J 

In the general case p > 3 a dynamical instability ap- 
pears at temperatures well above any relevant thermo- 
dynamic temperature. In particular, for p — 3, nu- 
merical integration of the dynamical equations as well 
as the use of series expansions (see below) yields a dy- 
namical transition at T)(p = 3) = 1.04 ± 02 in agree- 
ment with the inequalities (|l5|). Note that To is much 
higher than T d {p = 3) = 0.6125 or T s {p = 3) = 0.5. 



The relaxation time r re i ax associated to the decay of the 
distance D(t) to zero diverges according to a power law 
r relax ~ (T - r )-T with 7 ~ 1.1 ± 0.1. 

It is important to note that T is not related to any 
thermodynamic singularity. In the large p limit the 
inequality ([If]) yields To — » which gives a tem- 

perature much above the TAP temperature where an 
exponentially large number of states start to appear 
{Ttap —* y/\og{p)) ■ Consequently, the origin of the dam- 
age spreading transition is purely dynamical and not re- 
lated to any thermodynamic singularity or even to the 
existence of an exponentially large number of metastable 
states in the system. 

Now we discuss the behavior of the asymptotic dis- 
tance below T . In principle a new transition at T d 
(which we identify as Tl) is expected for the dynami- 
cal behavior of Q d {t) because the correlation C develops 
the mode-coupling instability. We will see that for p > 2 
the transition temperature T\ is absent. It is very diffi- 
cult to get analytical results below T d . A possible way 
to investigate the asymptotic long-time limit of Q d in 
the low-temperature regime (i.e below To) is to numer- 
ically integrate the set of dynamical equations (|7])-([ll]). 
Unfortunately, the CPU time and the memory needed 
to numerically integrate them grows very fast with the 
maximum time t (approximately like t 2 ). On the other 
hand, Q d (t) displays in several cases a non monotonic 
behavior as a function of time. Consequently, it is very 
difficult to extrapolate the numerical data to the infinite 
time limit. 

An alternative method was recently proposed in |Tl] ] 
were the series expansion for correlation and response 
functions was used to investigate the asymptotic long 
time limit of quantities such as the internal energy. Here 
we follow Jll| but extend their method within the con- 
strained formalism to include the series expansions for 
the correlation function Q(t, s). To this end, we decom- 
pose in Taylor series the correlation, the response as well 
as the overlap, 



(16) 
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where c k o = r k0 = S k0 and Q d (t) = J2kLo 9feo^ fe - It is 
possible to write in this case the recurrence relations be- 
tween the different coefficients c k i , r k i , qu , fJ-k ■ The time 
necessary to calculate the first coefficients of the series is 
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not very large and takes a few hours in a work station to 
reach the first 80 terms of the series Q The radius of con- 
vergence of these series is quite small. To enlarge their 
radius of convergence we have used Pade approximants 
to get an estimate of the asymptotic value of the distance 
Qd(oo). While the method works very well in case of the 
asymptotic value of the energy jllj (which depends on the 
Lagrange multiplier /i via the relation ju, = T — pE(t)) it 
is less effective for the asymptotic distance. The reason 
is that, while ji{t) is a monotonous increasing function of 
time, Qd(t) is not. In fact, below To the overlap Qd(t) 
has a maximum as a function of time for some values 
of T and the initial condition Qd(0). Consequently, the 
complex function Qd{z) turns out to have zeros close to 
the real axis and the radius of convergence of the Pades is 
smaller. But still it is possible to obtain some estimates 
for the asymptotic distance. Numerical integrations of 
the dynamical equations have been used to check that 
our extrapolations in the infinite time limit are correct. 




FIG. 1. Q d (t) with Q d (0) = as a function of 
time for different temperatures. From top to bottom 
T — 0.9,0.7,0.5,0.3,0.1. The continuous lines are the nu- 
merical integrations with time step At = 0.01 and the dashed 
lines are the reconstructed functions obtained from the Pade 
analysis. 

Some of our results are shown in figures 1 and 2 for 
case p — 3. We have studied three different initial 
conditions: a) anti-correlated random initial conditions 
with Qd{0) = — 1, b) uncorrelated random initial con- 



ditions with Qd(0) — and c) partially correlated ran- 
dom initial conditions with Qd{0) — 0.5. Case a) was 
analyzed using diagonal and the first off-diagonal Pade 
approximants assuming an asymptotic power law decay 
Qd{t) = Qd(oo)+ylt _7 . Cases b) and c) turned out to be 
more difficult to analyze due to the small radius of con- 
vergence of the series as well as to the presence of poles 
in the Pades. 

The behavior of Qd(t) for case b) is shown in figure 
1 for different temperatures. The continuous line cor- 
responds to the numerical integration of the dynamical 
equations while the dashed lines are the reconstructed 
functions Qd(t) obtained from the Pade analysis. Note 
the presence of a maximum in Qd(t) for several different 
temperatures. This feature is a consequence of the non- 
linear character (in Q) of equations (f||ll]) for p > 3 and 
is absent for p = 2 |l0| . 

Figure 2 shows the asymptotic distance for cases 
a),b) and c) as a function of the temperature. We find 
that the asymptotic distance is independent of the initial 
correlation. This is an interesting result since one would 
expect (at least below T<j) a dependence on the initial 
condition. We have to note that the dependence on the 
initial conditions is expected for models with the sym- 
metry a — » — er and with a simple free energy landscape. 
If the initial distance -Do is not zero there is always a 
finite probability (depending on Do) that the initial con- 
figurations start in different (symmetry related) ergodic 
components. Consequently, the asymptotic distance de- 
pends on the initial value of Dq. Here, such behavior is 
only found for p — 2 (l0| . 




*This is true for case p = 3 while for larger values of p the 
computational effort is larger 



FIG. 2. Asymptotic distance D x for p = 3 obtained from 
the Pade analysis of the series expansions for different initial 
conditions Do = 1 (circles), Do = 0.5 (triangles), Do = 0.25 
(stars). Typical error bars are shown for the last case. 
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In summary, we have studied the spreading of dam- 
age in the off-equilibrium mode-coupling equations. We 
have explicitly shown the existence of a damage spread- 
ing transition To and also found lower and upper bounds 
for its value. This transition takes place at very high 
temperatures. On the other hand, this transition is com- 
pletely unrelated to the existence of metastable states in 
the system. In fact, we have observed that this tran- 
sition is quite stable to the inclusion of any degree of 
asymmetry. Indeed in the fully asymmetric case we find 
that T ~ 0.71 for p = 3. Consequently, the damage 
spreading transition persists in the absence of the spin- 
glass phase or even in the absence of metastable states. 
This result corroborates some results already found for 
other disordered spin-glass models [jl2|,^3| . Interestingly, 
equations show that, only for p > 2, the coupling 

between the damage Qd(t) and the two-time correlation 
function Q(t, s) is non linear. This non-linear coupling 
is crucial for the appearance of the damage transition To 
which is well above Tj. We have also shown that below Tq 
the asymptotic distance is independent of the initial dis- 
tance. This was unexpected since such a dependence has 
been usually found in numerical investigations of several 
spin-glass models 14 1. It would be interesting to under- 
stand whether this result is a direct consequence of the 
first order nature of the glass transition. This and other 
issues, such as the scaling behavior of the overlap Q(t, s) 
in the aging regime and the existence of this transition 
in non-disordered glass forming liquids arc left for future 
investigations. 
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